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Abstract: We investigate signatures of the minimal supersymmetric inverse seesaw model 
at the large hadron collider (LHC) with three isolated leptons and large missing energy 
(3^ + $t ° r 21 + It + fx, with I = e, fj,) in the final state. This signal has its origin in the 
decay of chargino-neutralino (x\ X2) pair, produced in pp collisions. The two body decays 
of the lighter chargino into a charged lepton and a singlet sneutrino has a characteristic 
decay pattern, which is correlated with the observed large atmospheric neutrino mixing 
angle. This correlation is potentially observable at the LHC by looking at the ratios of 
cross sections of the trilepton + $t channels in certain flavour specific modes. We show 
that even after considering possible leading standard model backgrounds these final states 
can lead to reasonable discovery significance at the LHC with both 7 TeV and 14 TeV 
center-of-mass energy. 



Keywords: Beyond Standard Model, Supersymmetric Standard Model, Neutrino] 
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1. Introduction 

Supersymmetry is one of the most well-motivated theories for explaining new physics be- 
yond the standard model (SM) of particle physics. With the initiation of the large hadron 
collider (LHC) experiment at CERN, discovery of weak scale supersymmetric particles is 
highly envisaged. Moreover, the stable, neutral, and weakly interacting lightest supersym- 
metric particle (LSP) in an R-parity conserving theory can be an ideal candidate for a 
non-baryonic dark matter for explaining the data from the WMAP satellite and large scale 
structure formation. 

On the other hand, experimental evidences of neutrino oscillation have firmly estab- 
lished the necessity of massive neutrinos and the associated non-trivial mixing in the lep- 
tonic sector (see references [1-3] for the latest update on the three flavour global neutrino 
data). This is also a definite and very important indication of new physics because non-zero 
neutrino masses and mixing are not included in the SM. 

Seesaw mechanism [4-10] turns out to be the most simple and natural [11,12] way to 
accommodate tiny neutrino masses, m Ui (Yl m ^i < 0-58 eV [13]) both for supersymmetric 
and non-supersymmetric theories. However, neutrino mass models with canonical seesaw 
mechanism cannot be tested directly since the associated heavy states (to ensure smallness 
in neutrino masses) are usually ~ 0(1O 15 GeV) and thus well beyond the reach of any 
ongoing collider experiment. 

The inverse seesaw [14-16] is one of the viable alternatives with a seesaw mechanism 
operational at the TeV scale, which is well within the reach of the LHC. By proposal inverse 
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seesaw mechanism relies on a new mass parameter, which breaks lepton number explicitly 
by two units, and at the same time small enough to ensure small neutrino masses without 
the requirement of any super heavy states. The smallness of this new lepton number 
breaking scale can arise either spontaneously [16] or dynamically [17]. 

In the supersymmetric inverse seesaw model, i?-parity [18-21] is conserved as a con- 
sequence of lepton number violation by two units. Thus the LSP is stable for this class of 
models and can be a viable candidate for the cold dark matter of the Universe. 

In this paper, we study signatures of the minimal supersymmetric version of the in- 
verse seesaw model at the LHC. This model can accommodate three flavour global neutrino 
data [1,2] with a pair of SM gauge singlet superfields, which carry lepton number. These 
singlet scalars (having mixing with the doublet sneutrinos) can be thermal cold dark matter 
candidates, because the neutrino Yukawa coupling in this model is ~ O(10 _1 ). In a super- 
symmetric inverse seesaw model just one pair of singlets can explain the observed neutrino 
experimental data [22]. On the other hand, a non-supersymmetric inverse seesaw model 
requires at least two pairs of singlets [23]. Apart from the neutrino masses and mixing, 
the supersymmetric inverse seesaw model has also been analysed earlier in the context of 
lepton flavour violation [22,24-27], leptogenesis [28-30], dark matter [31-33] etc. 

In the case when one of the singlet sneutrinos is the LSP, the phenomenology at the 
LHC can be very interesting. This is because in this model the neutrino Yukawa coupling 
is large and can lead to lepton flavor violating (LFV) coupling of the sneutrinos with a 
charged lepton and the chargino. These couplings are related to the observed neutrino 
mixing angles and hence by studying the collider signatures of this model it is possible to 
study the relation between neutrino physics and the physics at the high energy colliders. 

Supersymmetric particle searches from the 1 fb _1 data, collected by ATLAS and CMS 
for pp collision at center-of-mass energy, ^/s = 7 TeV, has found no significant signal over 
the expected SM background. In the context of the constrained minimal supersymmetric 
standard model (CMSSM), searches by ATLAS exclude squarks and gluinos with masses 
below 950 GeV [34] at 95% C.L. for some particular choice of other parameters. The results 
from CMS extend the mass limit to 1.1 TeV [35-37]. However, the third generation squarks 
can still be somewhat lighter, particularly in the context of a more general MSSM scenario. 
This is the reason, in this work, we choose to work with a spectrum where the squarks of 
the first two generations and the gluinos are very heavy (~ 1 TeV) and the electroweak 
sector is relatively light so that the lighter chargino and neutralinos can be produced at 
the LHC. If in the production or in the decay chain the lighter chargino (xf ) appears 
then it can have a decay into a charged lepton (I) via xf — > P 1 + N, where N represents 
the singlet sneutrino LSP. The ratios of the decay branching ratios into different charged 
lepton flavors can be shown to correlate with the neutrino mixing angles [22]. Our aim in 
this paper is to look at these correlations by studying the trilepton + fx signature from 
the associated production of the lighter chargino (xf) and the second lightest neutralino 
(X§) at the LHC. 

Similar correlations also appear in the decay of the LSP in the model of bilinear Im- 
parity violation [38-45], spontaneous R-parity violation [46] and in //z/SSM [47-49]. In 
these models correlations with neutrino mixing angles have been studied in various context 



- 2 - 



in the case of a neutralino LSP decays as well as for other LSPs including the chargino [44]. 
The final states discussed in these cases generally include multi — leptons + jets + Tfs? 
along with the presence of displaced vertices originating from the long-lived LSP. In the 
present paper we have studied this correlation in the decay of the NLSP chargino in minimal 
supersymmetric inverse seesaw model (MSISM) through the cleaner trilepton + $t final 
state in the absence of any displaced vertex. 

It is known that at least two non-vanishing neutrino masses are essential [50-52] to 
account for oscillation data [1,2]. In the minimal supersymmetric inverse seesaw model 
with one pair of singlets only one neutrino mass is generated at the tree level, while another 
neutrino mass is originating from the sneutrino-anti-sneutrino loop [53,54]. This feature 
is analogous to the .R-parity violating (see review [55] and references therein) models of 
light neutrino mass generation with bilinear terms [54, 56-64]. Issues of neutrino mass 
generation in MSISM have been addressed in ref. [22]. Similar issues have been discussed 
in refs. [24,31], but with three generations of singlets. A model of neutrino mass generation 
where the origin of neutrino mass is radiative and suppressed by inverse seesaw scale has 
been advocated in ref. [65] . 

In MSISM, the atmospheric neutrino mixing angle (#23) correlates not only with the 
ratio of the branching ratios of the lighter chargino (x±) decay modes but also with ratio 
of the branching ratios of lepton flavour violating decays, r/fi — > £ + 7 [22]. On the other 
hand, trilepton signals (3/, with or without tau lepton(s)) have been extensively studied for 
a long time as an important probe for supersymmetric models [66-77] (see also references 
[3,4,6,7] of ref. [68]). Besides, a hadronically quiet event like this always has the favour of 
reducible backgrounds. Moreover, multi- lepton signals have already been considered as an 
important probe for seesaw models [78-80]. Being motivated by these features together 
with the novel correlations mentioned earlier, we aim to perform a detailed analysis of 
trilepton (3£, 2£ + It) + $t signals for the MSISM taking into account possible SM 
backgrounds. 

As mentioned earlier, we search for the trilepton signatures, arising from the decay 
of x^Xi 1 pah- In our chosen parameter points the associated production of the lightest 
chargino with the next-to-lightest neutralino, pp — > X^Xi 1 + X can occur with a detectable 
rate at the LHC. In addition, the lighter chargino and the second lightest neutralino decays 
via two body leptonic modes with large branching ratios. The final state signal will produce 
three charged lepton and missing energy signature (3£ + $t or 2£ + 1t + $t), because of the 
presence of the stable singlet sneutrino LSP. It is interesting to note that in our analysis the 
lightest neutralino decays into a singlet scalar LSP and a light neutrino. Both of these 
decay products escape detection and thus x? can be thought of as a virtual LSP, which also 
yields the missing energy signature at an accelerator experiment similar to that by an LSP. 
We investigate three body final states like 3£ (£ = e,n) and also 2£ + r-jet. Final states 
with more than one tau lepton have been dropped for small r detection efficiency [81]. 
In the course of present analysis we choose to work with non-universal gaugino masses 
but maintain M2 > Mi, where M\(M2) are the soft masses for U(1)(SU(2)) gaugino(s). 
It has also been assumed that \x > M2, where \i is the coefficient of the only bilinear 
term in the superpotential of the minimal supersymmetric standard model (MSSM). With 
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such a choice, xitXi are essentially gaugino like. Moreover, since the first two generation 
squark masses are heavy, the process pp —> X2X1 receives prime contributions from W^- 
boson mediated processes. Three of our benchmark points (BP1, BP3 and BP4 as defined 
later) are chosen with this criteria. We have, however, also considered the situation when 
M2 > n for another benchmark point (BP2). However, in a situation like this, X^iXi are 
higgsino like and consequently yield a smaller cross section for the process pp — > X2X1 • We 
show later that in a scenario like this the trileptonic final state possesses lower significance 
compared to the \i > M2 scenario. We will discuss this issue in more details later in section 
|5[ Having heavy squarks is also useful for suppressing flavour violating processes in the 
quark sector. 

The ratio of the branching ratios for Xi decaying into fj, and r channel in MSISM 
shows sharp correlation with tan 2 623 [22]. It is clear that one of the three charged leptons 
appearing either in 31 + fx or 21 + r + $t final states must have its origin in xf decay. 
Using this idea we find that the ratio ^(r^+f^aj > w ith ^ = e, /i shows nice correlation with 
tan 2 #23 even after the application of different kinematical cuts to reduce SM backgrounds. 
Definitely, that \i and r are coming from lightest chargino decay. Existence of this final 
state correlation with neutrino mixing angle along with a large amount of fix provides a 
distinct signature for the MSISM. 

The present paper is organized as follows. We start with a brief description of the 
underlying model in section || and discuss the generation of small neutrino masses as well 
as masses and mixing in the sneutrino sector. In section ^ we discuss the decay modes of 
the lighter chargino and the lighter neutralinos (XijX^) m our model and present a set of 
four benchmark points (BPs) studied in this work. In addition we discuss how these decay 
modes can lead to the trilepton final states at the LHC. Section |I] discusses the details of 
the signal event generation and the background analysis. We present our results in section 
H| and finally, our conclusions are provided in section ||. 



2. Model 

In the MSISM the particle content of the MSSM is extended by a pair of SM singlet 
fields, v c and S having lepton numbers —1 and +1, respectively. The model superpotential 
following refs. [22,31] is written as 

W = Wmssm + EabKL^Ht + M R v c S + ^fi s SS, (2.1) 



where Wmssm is the MSSM superpotential. In eq. (2.1) LjS are the SU(2)x doublet 
lepton superfields and H u represents a up-type Higgs superfield. v c represents a right 
handed neutrino superfield whereas S is another SM gauge-singlet superfield, but with 
non-zero lepton number. In eq. (|2.l| ) coefficient of the lepton number violating term is 
given by n$. In the limit us 0, MSISM superpotential (see eq. (|2.1|) ) restores lepton 
number conservation, which is consistent with the t'Hooft naturalness criteria [82]. 
The corresponding soft supersymmetry breaking Lagrangian is given by 

r — /-MSSM 1 _2 ~cf~c 1 _2 of q 
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e ah ^ h L\v c H h u + B Mr u c S + -B^SS + h.c. 



(2.2) 



with 



-£^| SM representing soft terms of the MSSM. Just like the coefficient MS appearing 
in the superpotential (eq. (|2.1|)), the 5 Ms parameter in the soft terms (eq. (p.2[)) violates 
lepton number by two units. The MSISM thus includes two lepton number violating 
parameters. Both of these will contribute to the Majorana neutrino mass matrix [22]. 

Tree level neutrino mass matrix in the MSISM is a 5 x 5 matrix in the basis (V/, v c , S), 
where, I = e, /U, r. This is given by 



/ m Dl \ 

m D2 

m D3 

m Dl m D2 m D3 Mr 

V Mr (1 8 J 



(2.3) 



where mo i = h\,v u (i = 1,2,3) are the three light neutrino Dirac masses with v u as the 
vacuum expectation value of the up-type Higgs field, {H®). The quantities h\, are the 
neutrino Yukawa couplings. The structure of the matrix shown in eq. fl2.3| ) can be readily 
understood from eq. ( |2.ip by looking at the mixing between different doublet and singlet 
neutral fermions. The effective 3x3 mass matrix for the light neutrinos can be obtained 
as 

Ms 



m Di m D] 



(2.4) 

becomes only 



tee ^ ( M R + Y, m D k ) 

In the seesaw approximation (m^ << Mr), the denominator of eq. (^ 
Mj|. The five mass eigenstates of eq. ( fE3| ) are denoted by rii, % = 1, . . . , 5, out of which 
^1,2,3 are nothing but three light neutrinos. 

Structure of eq. (|2.4| ) tells us that only one neutrino is massive at the tree level with 
mass: 

»'/,, (2-5) 



Ms 



The smallness of the neutrino mass is ascribed to the smallness of the fi s parameter, rather 
than the largeness of the Majorana-type mass, Mr, as required for the standard seesaw 
mechanism [83] . With the choice of normal hierarchy in the light neutrino masses this tree 
level mass will attribute to the atmospheric scale ~ 10~ n GeV. In the regime of seesaw 
approximation, for a typical Dirac mass m£> i ~ 10 2 GeV (assuming neutrino Yukawa 
couplings, hi ~ 10 _1 ) and TeV scale Mr, the value of the parameter fj, s comes out to 
be ~ 10~ 9 GeV. On the contrary, when rri£) i ~ Mr (see eq. (2.4)) then the atmospheric 
neutrino scale (~ 10~ n GeV) is determined by fi s only. 



The neutrino mass matrix shown in eq. (2.4) is diagonalizable using a 3 x 3 unitary 
matrix as 

U trT M v tree U tr = diag{0, 0, m U3 ). (2.6) 

The matrix U tr contains information about the tree level neutrino mixing angles. In order 
to satisfy the neutrino experimental data [1-3] one requires a second non-zero neutrino 
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mass eigenvalue which can arise by including the one-loop corrections in the neutrino mass 
matrix [22]. 

In this model the doublet and singlet sneutrinos mix after the electroweak symmetry 
breaking. Thus the sneutrino mass squared matrix, M~ is now a 10 x 10 matrix for MSISM, 
and assuming CP conservation this matrix can be decomposed into two 5x5 block matrices 
corresponding to CP-even and CP-odd sneutrino fields. The sneutrino mass term in the 
Lagrangian, then looks like, 



c, = -^ r a l ) + »; 2 r, . ( 2 - 7 ) 




where, (f> R = (v R , v cR , S R ), cj) 1 = (v- , u cI , S 1 ). The two mass squared matrices M± are 
given by [31] 



( (Af| + \M\ cos 2/3 + mjj.)Sij H A {vu - H m D] cot /3) m D] M R \ 

±( A h v v u - A» m A cot p) mlc + M R + J2 m 2 Dk fi S M R ± B Mr 

\ m Di M R hsM r ±B Mr ml + nl + Ml + B^J 

(2.8) 

where M~ denote soft supersymmetry breaking mass squared terms for SU(2)l doublet 
sleptons and Mz is the Z-boson mass. The ratio of the two Higgs VEVs is defined as 
tan/3 = 2st, where is the vacuum expectation value of the down- type Higgs field H ( 



The real symmetric mass matrix of eq. (2.7) can be diagonalized by a 10x10 orthogonal 
matrix as follows 

G M~ G T = diag(m 2 ^ , . . . , m%J, (2.9) 

< . . . < m~ . Di 

rately by 



"JVi'- ' ' ' "°Nio' 

with m ~ < . . . < m ~ . Diagonalizing the CP-even and CP-odd mass matrices Ml sepa- 



G± M| = diag(m~ ), i = l,...,5, (2.10) 

1 v izb 

where A^j + and iVj_ denote the z-th CP-even and CP-odd sneutrino mass eigenstates, 
respectively, leads to a different parameterization which can be used in some cases. In this 
notation, for the set of chosen parameters (shown later) 7Vi+ = A^2 and iVi_ = N\ and so 



on (see, eq. (2.9)). 



3. Decays of chargino and neutralino 

In this section we discuss the decays of charginos to charged leptons and singlet sneutrinos 
as well as the decays of the lighter neutralinos. We shall also show how these decays can 
lead to the final states, that we have proposed to study in this paper. Our choices of the 
four benchmark points for a detailed collider study will also be presented here. 

3.1 Chargino decay 

For the discussion of chargino decays we shall concentrate on a part of the parameter space 
where one of the singlet scalars of MSISM is the LSP. Hence this scalar singlet will appear 
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at the end of the supersymmetric cascade decay chains. For the present discussion let us 
assume that the dominant decay mode of the lighter chargino is in the two body mode 



Xi^N a + lt, a = 1,2, k = e,n,T, (3.1) 

with Ni being the CP conjugated state to N2- The relevant piece of the Lagrangian for 
the calculation of this decay width is 



Ltz-v = X-{Cl ja PL + C* ja P R )kN a + h.c. , (3.2) 



where 



Cfja — ~^/|b V jl(Gaj - iG a ,i+5) - KV* 2 (G a A ~ iG a g)}, 



c ija = ^=^iUj2(G ai -iG a ,i+b)- (3.3) 

The Y^s are the charged lepton Yukawa couplings and U, V are two unitary 2x2 chargino 
mixing matrices such that U*m 2x 2V _1 = diag(m~±,m~±), where m~±,m~± are the two 
physical chargino masses. The 2x2 mass matrix rri2 X 2 in the charged gaugino-higgsino 
basis ip + = —iX^Hu, = * s gi yen by 



"1-2x2 = • (3.4) 




Here g is the SU (2)l gauge coupling. 

The corresponding decay widths are given as 

r(xf -> Na + If) = % mn / a (I C ka I' + I C« a ft. (3.5) 

xt 

The members of CP conjugated pair of sneutrinos being nearly mass degenerate (m^ ~ 
m^) they are unlikely to be distinguished experimentally. Hence we sum over the CP-even 
and CP-odd sneutrino states of the CP conjugated pair. Thus 

2 

r(x± -+ iV 1+2 + If) = T (xt -+N a + If). (3.6) 

a=l 

One can adjust the parameters fig and B^ s in such a way that the tree-level neutrino mass 
matrix contribution determines the atmospheric mass scale, while the one-loop corrections 
control the solar mass scale [22] . In such a situation it can be shown that in order to have 
small reactor neutrino mixing angle and maximal atmospheric neutrino mixing angle, the 
parameter mr> 1 has to be considerably smaller than other two Dirac masses and simul- 
taneously, iriD 2 ~ iriD 3 . The solar neutrino mixing angle can be kept large by keeping 
the parameters Si = A l h ^v u — fim£> i cot/3 to be of the same order for all the three flavors, 
i = e, fi, t. In this case, one can show that the decay width of the lighter chargino, 
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T(xf — > N\ + 2 + li ) correlates with the corresponding parameter nip.. The atmospheric 



9 _D 

neutrino mixing angle at the same time also behaves as tan #23 ~ — r 2 "- Hence, one would 

m D 3 




Figure 1: Correlation plot for the ratio of the branching ratios Br< -^ ~^ I ^ 1 + 2+fl } with —S 2 - 



3.2 Neutralino decay 

In our chosen benchmark points (defined later in this section) the lightest neutralino is the 
next-to-lightest super symmetric particle (NLSP) and, decays dominantly through the two 
body decay channels Xi ~^ u k + ^1,2, h = e,fi, r. The relevant interaction term of the 
Lagrangian is: 

L U x°u = xM L mjb P L + A* jb P R )is m N b + h.c, (3.7) 



where 



Amjb = § ( N j2 -tan6 w 'N* 1 )(G b i - iG fe ( i+5 ))£/* , 



Amjb — -^K U im^j4:(G b 4 - iG b9 ). (3.8) 

Here g is the SU(2)l gauge coupling, 9w is the weak mixing angle, and N is the unitary 
4x4 neutralino mixing matrix. Although the second lightest neutralino (jq) decays mostly 
through the standard MSSM two-body charged lepton-slepton channel (X2 ~ * h + 
some of its branching fraction goes into the decay channels arising from the coupling given 
in eq. ( |3.7[) . Here, we have neglected the charged lepton flavor violating decay of x%- The 
decay width of a neutralino decaying into neutrino-sneutrino two-body mode is given as 

_ (m 2 — m 2 ~ ) 2 

T&j N b + u m ) = |— ^ (| A L mjb | 2 + I A* jb | 2 ). (3.9) 
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3.3 Trilepton signal and the benchmark points 

In order to illustrate the trilepton signal we simulate X2X1 production followed by their 
two-body decays to produce 3£ + $t or 21 + r — jet + $r final states, where £ = e, /z. 

M ±/r± 



.^1,2 




Figure 2: Feynman diagram for the process pp — > 3£ + Ifir or 2£ + t — jet + Ifix- 

As discussed above the production process and the decay cascades leading to these 
final states are as follows 



PP - 

xt 
xl- 

7± 



X°2 + Xf, 

£ ± +£ T , 
ui + N h2 . 



(3.10) 



The Feynman diagram for the above mentioned final states is shown in figure ^. In the 
presence of heavy squarks (~ 1 TeV), this is the leading process for the chosen signal. 

Because of the presence of the massive singlet sneutrino LSPs, N± t 2, (quasi-degenerate 
in masses), we have, for this model, a large amount of missing energy in the final states. In 
order to have an appreciable signal rate one must have significant production cross section 
of X2~ xf P a h an d large branching ratios for the above-mentioned decays. To achieve these 
we have chosen four benchmark points (BPs) in the parameter space where the detailed 
collider simulation has been performed. We scanned the whole parameter space to check 
for charged lepton flavour violating (LFV) decay widths and we found points both above 
and below the experimental limits in different region of the parameter space. In all of the 
benchmark points, constraints from LFV decays [84] are satisfied as well as the atmospheric 
neutrino mixing is near maximal. The input parameters for different benchmark points are 
given in table |]. The choices of the parameters m? D _ will be shown later. 

The mass splittings between the second lightest neutralino, the charged sleptons and 
the lightest neutralino are maintained in a way, that the second lightest neutralino decays 
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BP1 


BP2 


BP3 


BP4 


tanp 


5 


10 


20 


10 


H (GeV) 


330 


240 


280 


350 


Tl AT / S"-< T 7"\ 

Mi (GeV) 


170 


195 


160 


240 


M2 (GeV) 


220 


340 


240 


290 


M 3 (GeV) 


1100 


1100 


1100 


1100 


Mr (GeV) 


145 


IfaO 


14U 


150 


fi s x 10 (GeV) 


1 .80 


/ .81 


1 . ID 


1 . lb 


to 2 (n e y^) 

1 1 VyC \ VJC V J 


2500 


3025 


2500 


3025 


m| x 10" 4 (GeV 2 ) 


4.0 


4.0 


4.0 


4.8 


£ A/fl (GeV 2 ) 


2500 


2500 


3500 


2500 


S MS (GeV 2 ) 


10 


10 


10 


10 


M| x 10- 5 (GeV 2 ) 


5.63 


5.63 


5.63 


5.63 


M| x 10- 4 (GeV 2 ) 


2.99 


3.69 


2.59 


5.86 


M4 x 10- 4 (GeV 2 ) 

e 22 


2.99 


3.69 


2.59 


5.86 


M 2 x 1(T 4 (GeV 2 ) 


3.53 


7.08 


7.90 


8.18 



Table 1: Values of the relevant input parameters for different benchmark points. The quantities 
M| c represent soft squared masses for the right-handed charged sleptons. Other parameters are 
defined in the text. 

only through charged lepton-slepton two body modes and the charged sleptons further 
decay into the lightest neutralino and charged lepton states. With these considerations we 
generated the sparticle spectrum using SuSpect (version 2.41) [85]. Masses of the neutrino 
and sneutrino states are computed using a self developed code in FORTRAN. Relevant mass 
spectra for these benchmark points are shown in table ^. 

The choice of model parameters for different benchmark points are chosen to yield 
statistically significant final states. As an illustrative example, production cross sections 
for the x^x^ 1 pair with 7 TeV center of mass energy at LHC are in the range of 200 — 300 fb 
for the first and third benchmark points. For the fourth benchmark point with relatively 
heavy x^xf P arr ( see table |2|) the production cross section is reduced by a factor of 4(3) 
compared to the first(third) benchmark point. On the contrary, a higgsino like x^xf P an " 
(BP2) yield a similar production cross-section like BP4, in spite of having a lighter X2X1 
pair. Thus, the region of parameter space with higgsino like x^Xi P a i r is unfavorable for 
this analysis. 

Note that hadronically quiet trilepton signal (3£ + $t) will get very little contribution 
from squark-squark, squark-gluino and gluino-gluino pair production. On the other hand, 
when we have 2£ + r-jet + Tfs? signal, then one should consider all other sources of dilepton 
+ 1-jet + fx events where one jet can be faked as a r-jet. For example, one can have a 
jet out of a squark decay (q — > q' + %\ ) from one side of the cascade. However, since in 
this model the squarks are much heavier (~ 1 TeV) and after incorporating the probability 
of any jet faking as a r-jet, the event rate comes out to be negligibly small compared to 
the one generated from chargino-neutralino production. Hence the main contribution to 
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2£ + r-jet + $t signal comes from X2X1 production only. 



4. Event generation and 
background analysis 

On the basis of the discussion presented 
in the previous section, let us now provide 
a detailed description of event generation 
and subsequently, the background analy- 
sis. The decay widths corresponding to 
the two-body modes shown in eq. ( 3.10Q 



BP1 



BP2 



BP3 



BP4 



N 1+ 


153 


27 


169 


18 


147 


96 


159 


76 


N 2+ 


247 


38 


256 


45 


244 
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266 


53 
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27 


169 


18 


147 


96 
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76 


N 2 - 


247 


34 


256 


42 


244 


76 
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B4 


145 


46 


160 


52 


140 


42 


150 


45 


n 5 


145 


46 


160 


52 


140 
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eL,fiL 


751 


28 


751 


37 


751 


39 


751 
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00 
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166 


83 


245 


98 
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193 


11 


269 


48 


284 


12 


289 


17 




751 


30 


751 


40 


751 


53 


751 


43 


X? 


159 


51 


172 


90 


151 


80 


226 


36 




198 


11 


234 


59 


207 


67 


264 


16 


~± 

Xi 


192 


79 


215 


47 


203 


72 


255 


99 


~± 

x 2 


363 


91 


372 


89 


326 


77 


391 


60 



Table 2: Relevant mass spectra obtained for four 
benchmark points with Tii,2,3 = ^1,2, 3- 



have been used to modify the branching 
fractions of the charginos and neutrali- 
nos obtained from SuSpect. These in- 
put files are then fed to PYTHIA (version 
6.409) [86] for event generation and show- 
ering. Initial and final state radiation, de- 
cay, hadronization, fragmentation and jet 
formation are implemented following the 
standard procedures in PYTHIA. Factoriza- 
tion and renormalization scales are set at 
yfs (i.e fj,R = fiF = ), where is 

the parton level centre of mass energy. We have used the leading order CTEQ5L parton 
distribution functions [87,88] for the colliding protons. Some of the background events are 
generated using ALPGEN (version 2.14) [89] with default factorization and renormalization 
scales. The jets are constructed using cone algorithm in PYCELL. Only those jets are 
constructed which have pt > 20 GeV and | rj |< 2.5. To simulate detector effects we have 
taken into account smearing of jet energies by a Gaussian probability density function of 
width [90] cr(E)/Ej = (0.6/ \/E~\GeV}) + 0.03 where Ej is the unsmeared jet energy. 

In order to find three isolated leptons in the final states we impose following cuts and 
isolation criteria: 

I. Leptonic events are selected only if p T > 8 GeV and | rf \ < 2.4. 

II. Lepton-lepton separation AR(£,£) set to be > 0.2, where AR = \J (Arj) 2 + (Acj)) 2 . 

III. Lepton-jet separation AR(£,j) chosen to be > 0.5. 

IV. The sum of Et deposits of the hadrons which fall within a cone of AR < 0.2 around 
a lepton, must be less than 10 GeV. 

A pt cut of 10 GeV and 17 GeV [81] is applied on final state muons and electrons 
respectively, for the analysis at 7 TeV and 14 TeV center of mass energies at the LHC. The 
r-jets are counted withp-p > 20 GeV and | if |< 2.4. The r's are then counted according to 
the visible energy bins. A T-jet is treated as tagged or untagged according to the efficiency 
(e r ) of the most efficient algorithm given in [91]. In reference [91], r identification efficiency 
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obtained from actual collision data at 7 TeV center of mass energy has also been quoted. 
The efficiencies obtained from Monte-Carlo simulation and from the data agrees very well. 
However, for higher luminosity with 14 TeV center of mass energy, a lot of underlying 
events are expected to be there, which can perhaps bring down the detection efficiency. In 
this case also we have used the same efficiency as in 7 TeV case hoping the experimentalists 
can maintain the efficiency as we have now. Unlike r, detection efficiencies of e and [i are 
assumed to be 100%. 

We have analysed the SM backgrounds in some detail. The dominant background 
events arise from tt and WZ production at the LHC. Apart from these, contributions 
from ZZ, WW, Zbb, Wbb, Z + jets, Wt, tb, WWW, Wtt events have also been studied 
at the leading order. We also studied QCD di-jet events. But after putting the cuts to 
reduce backgrounds as mentioned below we found no trilepton events for 1 fb _1 integrated 
luminosity from these particular QCD events. We use ALPGEN for an estimation of Zbb, 
Wbb, Wt, tb, Z + jets, WWW, Wtt backgrounds. We generate these events at the parton 
level using ALPGEN and fed those partonic events to PYTHIA for showering, hadronization, 
fragmentation, decay, etc. The other events are generated and analysed using PYTHIA. It 
should be mentioned that the importance of these processes have already been emphasized 
in the literature [74,92]. 

The trilepton signal in our model arising out of chargino-neutralino production is 
accompanied by large missing transverse energy ($t) > because of a pair of singlet sneutrino 
LSPs and a neutrino. As an example, the fix spectrum of background events as well as 
the signal events (3/U + $t) for the first benchmark point (BP1) are shown in figure |||. 
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Figure 3: Ifir plot for signal events 3/j. + $t and summed up contribution coming from the 
WW, ZZ, Zbb, Z + 1 jet, Z + 2 jets background for LHC at 14 TeV center of mass energy with 1 
fb _1 of integrated luminosity. The $r bin size is chosen to be 10 GeV. 



These distributions are obtained without applying any cuts to reduce background 
events. It is evident from the plot in figure || that a strong I^t cut will affect the signal 
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cross-section very mildly, but it reduces significantly background events coming from some 
processes. Therefore, a cut Jfix> 25 GeV is applied for background rejection. For some other 
channels; tt, WZ, Wti, WWW the fx distributions do not peak before 25 GeV as shown 
in figure Hence, the above mentioned fix cut does not seriously affect these background 
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Figure 4: $ T plot for 3^ + $ T events obtained from tt, WZ, Wtt, WWW backgrounds for LHC 
at 14 TeV center of mass energy with 1 fb _1 of integrated luminosity. The fx bin size is chosen 
to be 10 GeV. 



events. To reduce these events we have further applied two more cuts. An invariant mass 
cut on the opposite sign dilepton pair, 80 GeV > M?~ v > 100 GeV removes backgrounds 
coming from Z-bosons. To manifest this idea we show invariant mass distribution in figure 
|5] constructed from opposite sign muon pairs for signal events and WZ background events. 

10° C , 
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Figure 5: M w plot for signal events and WZ background for LHC at 14 TeV center of mass 
energy. The M w bin size is chosen to be 5 GeV. 
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On the other hand, rejection of tagged 6-jet events significantly reduces backgrounds 
coming from tt events. A jet (with | r] |< 2.5) is reconstructed as a b-jet if the AR 
separation between the jet and the 6-quark (with px > 5 GeV) is less than 0.2. The b-jet 
identification efficiency is taken to be 50%. 

In order to perform the collider analysis we have randomly generated m 2 D _ and Sf 
within certain range: (Yli m D. l ) 1/2 G 10 [ ~ 4 ' 2 - 61 and (Ei^ 2 ) 1/4 G 10 [_4 ' 3] [22]. Moreover, we 
also consider i^2,fn 2 Di ) 1 / 2 > (J2 i 5 2 ) 1 / 4 , such that Dirac neutrino masses give the dominant 
contribution to the chargino decay [22]. Around each of the four benchmark points we 
select a set of six to seven points of these randomly generated parameters. These points 
will be useful for the correlation study discussed later in section ||. Remember that these 
parameters control the neutrino masses and the mixing angles and our choices of benchmark 
points are such that the atmospheric neutrino mass scale is determined by the tree level 
neutrino mass matrix contribution. Before showering in PYTHIA, as mentioned earlier, the 

ratio - Br ^i ~ > ^ 1+2+ ^ } shows a very nice sharp correlation when plotted against —w 2 - which 

Br( x f^N 1+2 +T±) v & m 2 D3 

is a measure of tan 2 #23- We have done the showering for four benchmark points introduced 
in table Q and table || to look for the ratio ^^Epj; S ggjga with £ = e, \i. Since one fi and 
one r in these final states always come from the decay of xti we would expect that this 
ratio will also go as ~ tan 2 #23- Hence, by measuring this ratio from the trilepton signals 
one can obtain information about the atmospheric neutrino mixing angle at the LHC. On 
the other hand, a precise measurement of the atmospheric neutrino mixing angle at the 
oscillation experiments can be used to predict the allowed range of the above ratio at the 
LHC. In the following section we give a quantitative estimate of this ratio for our choices 
of benchmark points (along with randomly selected values of m 2 D .) and show that for each 
of these points the various signal events included in the calculation of this ratio can be 
statistically significant. 



5. Results 

In order to study the correlation between the atmospheric neutrino mixing angle (#23) 
and the final states with trilepton + ^rat the LHC, we look at the ratio of cross sections 
ct(pp^t ee+fcr) ' ^ = e ' /U ' men ti° ne d in the introduction, in the denominator the r must 
always come from the decay of Xi because we are considering final states with only one 
r-jet and neglecting lepton flavor violating decays of X2 an< ^ For the same reason, in 
the numerator one fi must always also come from the decay of %\ ■ Hence, naively we would 
expect that this ratio of cross sections will also show nice correlation with the atmospheric 
neutrino mixing angle #23- 

After applying different cuts to reduce backgrounds and taking into account the r- 
tagging efficiency, we find that the ratio of trilepton signal cross section again shows a nice 
correlation with the atmospheric neutrino mixing angle tan 2 #23- However, in this case the 
numbers change from the ratio of branching ratios, discussed earlier and the straight lines 
obtained are steeper than the one shown in figure 0. This happens because in our simulation 
we take the detection efficiency of fi to be 100% as opposed to the r detection efficiency, 
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which is smaller [81]. Since the branching fractions of r events are in the denominator of 
the ratio, the numbers naturally go up. 

The cross-sections and the corresponding statistical significance ( tw^te with x = 
e, (j,, t) obtained from our simulation for LHC are shown in this section. Here S x is defined 
as the number of x ^ it signal events and B x is defined as the number of corresponding 
background events. In more simple form significance for the /x ^ ££ + $t channel is defined 

c< -\-S *i 

as mmm _ j n a similar fashion significance for the T^tt + $t channel can 

be obtained. 

We quote the results below for an integrated luminosity of 25 fb _1 for the LHC with 
7 TeV and 14 TeV center-of-mass energies. The results are obtained with the cuts men- 
tioned in section ^. Throughout this analysis we have used leading order cross sections for 
the signals as well as all the backgrounds at the LHC. However, if next-to- leading order 
(NLO) corrections are included the statistical significance will not change much. For ex- 
ample, if NLO corrections are included the signal cross section at 14 TeV LHC is expected 
to increase by 1.25 to 1.35 [93]. As discussed above, a large contribution to the background 
comes from the tt events. The NLO cross section for ti production at 14 TeV LHC is about 
800 pb [94, 95] which is about a factor of two larger than the leading order cross section 
that we have used in our analysis. Thus taking into account the NLO contribution of all 
the major background events along with the signal event, the significance S x /\/B x + S x 
estimated for our signal, will not change much and remains conservative in comparison to 
the uncertainties in the production cross sections. 

Values of the randomly generated parameters mfj., for four chosen benchmark points 
are presented in table |3[ For the numerical analysis we choose to vary m 2 D in the range 
of 10 -4 - 10 -2 GeV 2 , whereas m 2 n are varied within 10~ 2 to 10 2 GeV 2 . The 5} are also 
varied accordingly, but keeping the constraints > (S^) • The scale of m 2 D . 

has a strong influence on the decay processes xf — > Ni + 2 + ^ / T± and Xj — > + v m . In 
order to achieve a statistically significant trilepton final state originating from x^xf pair, 
we would like to have Br(xf — > + /i^/r^) to be large and Br(x% — * Nt, + z^ m ) to 

be small, simultaneously. However, in the limit m£>. ~ Mr ~ O (10 2 GeV), the neutrino 
Yukawa couplings h l u are ~ O (1). Then as can be seen from eqs. ( |3.5| ) and ( |3.9| ) both 
of these decay widths are large and consequently, yields a smaller branching ratio for 
X° - > & + P ■ We observe that in this case it is rather difficult to achieve a statistically 
significant final state particularly for the r it + fx mode. 

In the trilepton signals studied in 
this work, one lepton comes from the 
lighter chargino (xf ) decay and the other 
two same flavour opposite sign leptons 
come from the second lightest neutralino 
(X°) decay. Since the probability of get- Table 3: Randomly generated values of m%. corre- 
ting electrons from the chargino decay is spending to the four benchmark points as indicated 
suppressed compared to muons or taus, ^ n sec ^ 011 EI 

events with odd number of electrons (eee and should have smaller cross-sections 



m 2 Di (GeV 2 ) 


BP1 


BP2 


BP3 


BP4 


m 2 Di x 10 4 


2.12 


3.30 


1.58 


4.24 


m 2 D2 


71.80 


62.96 


68.50 


80.33 




62.45 


54.87 


66.63 


86.00 
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compared to others, which is clearly reflected in the signal cross-sections. This feature is 
intrinsically related with the small but non-zero reactor neutrino angle [96], which will be 
discussed again later. In table || and || chosen trilepton +f>T cross sections are shown along 
with the total standard model background cross section for the LHC at center-of-mass 
energy, yfs = 7 and 14 TeV, respectively. The corresponding statistical significance of the 
signals are shown respectively in table || and table [7|. We can see from table || that, at the 
LHC even with yfs = 7 TeV, the lowest signal significance for r ^ it + fx final state, that 
we have obtained, is greater than 3a for an integrated luminosity of 25 fb _1 . Hence, the 
trilepton +$t data for 25 ftv 1 integrated luminosity at 7 TeV LHC should be able to 
constrain the theoretical parameter space of this model. These numbers (significance) are 
much higher for LHC with yfs = 14 TeV and are shown in table [?]. It is once again evident 
from these tables that a higgsino like X2X1 P a ^ T (BP2) yields statistically less significant 
specific trilepton final state. In other words for such benchmark points, the significance of 
the final state trilepton signal is less promising. This situation is comparable to a heavy 
gaugino like X2X1 pah as represented by BP4. 



Tri-lepton 




signal 






Background 


events 
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<x(fb) 
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BP2 


BP3 


BP4 




eee 


0.37 


0.31 


0.50 


0.23 


8.73 


eefj, 


9.47 


5.37 


8.30 


3.12 


18.91 




1.08 


0.49 


1.26 


0.66 


21.15 


nnn 


24.13 


8.21 


18.85 


8.51 


23.84 


eer 


2.93 


2.14 


2.86 


1.40 


4.60 


HfiT 


7.17 


3.39 


6.99 


4.04 


13.18 



Table 4: Cross-section for different trilepton channels are shown here for four different benchmark 
points along with their total SM background contribution for LHC with yfs — 7 TeV. Corresponding 
input parameters and mass spectrum are given in table 0, table and table 0, respectively. . 



Tri-lepton 




Significance 






events 




















BP1 


BP2 


BP3 


BP4 


fxee + nfifi 


19.23 


9.05 


16.24 


7.89 


ree + Tfj,fj, 


9.57 


5.72 


9.37 


5.65 



Table 5: Statistical significance of the studied trilepton signals with integrated luminosity 25 fb 1 
at the LHC for yfs = 7 TeV for different benchmark points. 

We present the correlation plots, obtained with different randomly generated values of 
m 2 D , and 6f around each of the four benchmark points. These are shown in figure ^ and 
figure ^ for the LHC with yfs = 7 TeV and 14 TeV, respectively. We present these 
correlations with best fit lines. It can be seen from these figures that the 3a allowed value 
of tan 2 023 [1] from atmospheric neutrino oscillation experiments predict a value of the 
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/W 


65.42 


24.32 
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26.47 
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eer 
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6.48 
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4.76 


29.16 


/i/JT 


20.11 


9.80 


18.46 


12.48 


56.29 



Table 6: Cross-section for different trilepton channels are shown here for four different benchmark 
points along with their total SM background contribution for LHC with = 14 TeV. Correspond- 
ing input parameters and mass spectrum are given in table [l], table [| and table g respectively. 
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Significance 
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\/B x j2tt+Sx'E,tt 








BP1 


BP2 


BP3 


BP4 


/nee + fififi 


28.28 


13.55 


24.02 


12.61 


ree + t/i/j, 


13.53 


8.07 


12.97 


8.50 



Table 7: Statistical significance of the studied trilepton signals with integrated luminosity 25 fb 
at the LHC for y/s = 14 TeV for different benchmark points. 



ratio of cross sections ^ et+fcr) ' ^ = *° ^ e a PP rox i ma t e ly i n t ne range 1.0 — 6. 

These predictions can be verified at the LHC or the measured value of this ratio can give 
an alternative estimate of tan 2 #23- On the other hand, if this ratio comes out to be very 
much different from the ones predicted here then one can perhaps conclude that MSISM 
is not the correct model for explaining neutrino masses and mixing. 

Nevertheless, as we can see, from the correlation plots, that there is a different linear 
relationship for each different kind of benchmark points. In general then, from neutrino 
oscillation data we cannot give a unique prediction for the ratio of the cross-sections that 
can be verified at the LHC and help us in constraining the model parameters. In other 
words, measuring the cross section ratio at the LHC would not allow a prediction of 623 that 
could be tested against oscillation results. This means that we need other measurements 
at the LHC to allow such predictions. As an example, to distinguish among the four 
benchmark points we plot the ratio, m~± /m^^ 2 with the ratio of cross-sections of in and r 
channels which gives four separate parallel lines for the four benchmark points (figure 
One can see from figure [8| that the ratio { m ^± / m jy 12 ) increases as the slope of the straight 
lines in the correlation plot corresponding to different benchmark points decreases. This 
pattern can easily be understood. Increase in the mass ratio indicates greater splitting 
between the chargino and sneutrino masses. As the splitting increases, the leptons coming 
from this chargino decay become more energetic (eq.( 3.lti| )). This affects the r count in the 
final state more than the /x count as the detection efficiency for the taus increases with the 
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Figure 6: Correlation plot ( CT (pp_»^± ^ ii+fcr) vs * an 23 ' w ith £ = e,fx) obtained for LHC at 
y/s = 7 TeV. Three differently coloured vertical strips correspond to la, 2a and 3a allowed 
region for tan 2 6 l 23, respectively. The benchmark points as given in table [|, table || and table || are 
represented by coloured ■. Other points, represented by coloured •, are obtained with randomly 
generated m 2 D . and S\ values. 




increase of visible energy of r decay products [91]. Hence more r events are expected in the 
final state for those benchmark points which has greater lighter chargino - LSP mass ratio 
for a given set of m^'s. With the increase of r events the ratio of the cross-sections plotted 
in the correlation plots decreases and as a consequence gives smaller slope compared to 
the previous benchmark point. Now it is clearly understood that if we can determine the 
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Figure 8: m~±/m^ i 2 plotted against 

as the slopes of the straight lines in figure H 



corresponding to the four benchmark points decrease. 



for = 7 TeV. The mass ratios increase 



lighter chargino-LSP mass ratio, we can distinguish among the four benchmark points. 

Mass determination techniques in the context of LHC have been studied extensively. 
Transverse mass variable (mr 2 ) [97,98] is very useful for this purpose. rriT 2 has also been 
generalized for the cases where the parent and daughter particles in the two decay chains 
are not identical [99,100]. Moreover, final state with more than two invisible particles has 
also been addressed in ref. [101]. In our case, we observe the following: 



• The lightest neutralino Xi is also invisible, as mentioned earlier in the text. 

• One lepton is produced from one side of the cascade and remaining two leptons from 
the other side of the cascade (figure . 

• X\ and X2 are not mass degenerate but the difference is quite small in the context of 
mass measurement. 

So, we see that daughters of different masses are produced here from nearly identical 
parents. A mass determination technique similar to refs. [99-101] can be applied here too 
to determine the masses of the lighter chargino and the sneutrino LSP. However, a detailed 
analysis in this direction is beyond the scope of the present paper. Thus we see that 
measuring the mass ratio (m~± /m^ 2 ), along with the ratio of the trilepton cross-section, 
can help us pick the correct benchmark point and hence predict the correct value of #23 that 
could be tested against the oscillation results. On the other hand, a precise determination 
of tan 2 623 from oscillation experiments as well as a measurement of the cross section ratio 
at the LHC can give a unique prediction of the mass spectrum of the model, that can be 
verified by mass measurements at the LHC. 

In support of our explanation for obtaining different slopes, we present the following 
analysis. Since this difference among the four benchmark points appears because of taking 
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different r identification efficiencies for different energy range and for taking separate pt 
cuts for //'s and r's, we can remove this by the following strategy: 

• A pt cut of 20 GeV taken for both fi and r. 

• A uniform r identification efficiency of 50% applied over the whole energy range. 



12 




Figure 9: Correlation plot ( ct[pp^>t± fj ££+Pt) vs tan2 ^23) obtained for 7 TeV center of mass energy 
under the assumption of uniform t identification efficiency (e r ~ 0.5) and same cut for both fj, 
and r. 

We have presented the result in fig.||. This shows the correlation plot for 7 TeV center 
of mass energy under the above mentioned conditions. All the benchmark points now lie 
almost on one straight line. A few points still looks a little bit scattered because of the 
different isolation criteria used for /x and r. 

Finally, in table | we show the ratios ^p^se+fir) and ^^t+fo) for the four 
benchmark points. The smallness of these ratios is also a distinct feature of this model and 
arises due to the smallness of the neutrino reactor angle imposed by neutrino data [96]. 
In the usual MSSM scenario these ratios are expected to be much higher as there is no 
suppression of charginos decaying into electrons as we have in this model. 

A discussion of these specific trilepton signals remains incomplete without a note on 
the Tevatron analysis of the considered model. For the four chosen benchmark points we 
observed no points with significance > 3a for the r U+^T final state and simultaneously 
consistent with the atmospheric neutrino mixing at the 3a limit. This is a well expected 
result considering that the Tevatron center-of-mass energy is 1.96 TeV with 12 fb _1 of 
integrated luminosity [102]. For example, the statistical significance for mode for 

BP1 with m^.s given in table ^ is computed to be 1.64. 
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5.24 
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BP4 
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5.38 


1.89 



Table 8: The ratios °%^£+£$ and for the four benchmark points. 



6. Conclusion 

We consider the minimal super symmetric inverse seesaw model and study its characteristic 
signatures at the LHC. This model, with only one pair of singlet superflelds explains existing 
neutrino oscillation data. The model is rich from phenomenological point of view and can 
lead to potentially testable signatures at the hadron colliders. In this R-parity conserving 
model, one of the singlet sneutrino (with a small admixture of the doublet sneutrino) is 
the lightest supersymmetric particle (LSP) and as a result shows up in the collider as 
missing energy. Charginos can decay to charged leptons plus singlet sneutrino LSP. The 
decay patterns of the chargino are controlled by the same parameters which generate the 
neutrino mixing angles. 

In order to study this correlation of the chargino decays and the neutrino mixing angles, 
we look at specific trilepton + j$ T signatures at the LHC. We show that the ratios of cross 
sections of this studied trilepton + fl T final states in certain flavour specific channels 
(fiee + fir, fi^fi + fix, Tee + fir, Tfifi + $t ) nicely correlate with the atmospheric neutrino 
mixing angle. We explore different points in the parameter space to study this correlation. 
A measurement of these cross sections thus provide an interesting test of the minimal 
supersymmetric inverse seesaw model. The hard missing Et spectrum makes this trilepton 
final state statistically significant by reducing certain standard model background events 
significantly. We adhere to different cuts to reduce the backgrounds coming from some other 
channels. Motivated by the recent results from the ATLAS and the CMS experiments, we 
work in a scenario with heavy squarks and gluinos and a relatively light electroweak sector. 
The results of our analysis suggest that the theoretical parameter space of this model can 
be constrained by the data collected at the LHC with center-of-mass energy 7 TeV and 
for an integrated luminosity of 25 fb _1 . On the other hand, a measured value of this 
ratio at the LHC can give us an alternative estimate of tan 2 #23 and confirm (or rule out) 
this minimal supersymmetric inverse seesaw model as a possible explanation of neutrino 
masses and mixing. We also show, as a distinct feature of this model, the cross sections of 
pp —7- efifi + ftj< and pp — > eee + are suppressed compared to the total chosen trilepton 
+ fix cross section because of the restrictions on the neutrino reactor angle imposed by 
neutrino data. 
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